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Abstract 

By a weak deformation of the cylindrical symmetry of a potential vortex in a 
relativistic perfect isentropic fluid, we study the possible dynamics of the central line 
of this vortex. In "stiff " material the Nambu-Goto equations are obtained. 

1 Introduction 

These last years, the motion of the relativistic vortices has been in general studied in the 
context of a model based on nonlinear wave equations of a complex scalar field. These 
studies enlight the theoretical similarities between superfluid vortices and global cosmic 
strings §@,|,!§. 

In this work, we take a purely hydrodynamics point of view The purpose is to find the 
relativistic dynamics of the central line of a potential vortex in an isentropic perfect fluid 
when it deviates weakly of the stationarity and of the cylindrical symmetry. 

We suppose that the core of this vortex is thin and that its central line has a small 
deviation from the straight configuration; so that locally, in the neighbourhood of the core 
(which is empty), the flow must approach the cylindrical symmetry. Therefore, expanding 
the equations of motion in function of small parameters caracteristic of the curvature and 
the thickness of the core, the lower order which is the equation of the straight vortex, is 
automatically verified. The following order gives approximatly the dynamical equation of 
the central line. This method has been transposed from a study of the dynamics of a 
self-gravitating thin cosmic string ||. 

The plan of the paper is the following. In section 2, we shall review the essential 
ingredients of a perfect fluid and give the formulation of Lichnerowicz in which it is easy 
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to write the stationnary cylindrical solutions of an isentropic fluid in Minkowski spacetime 
. In section 3, by perturbating the cylindrical symmetry, we obtain the dynamical equation 
of the central line. The interesting result is that in stiff matter it is reduced to the Nambu- 
Goto dynamics. In section 4 we discuss the solutions and finally we end with a conclusion. 



2 Potential vortex in a perfect isentropic fluid 

Let us consider a perfect fluid ||. The tensor energy-momentum is 

= (p + p)unU„ + pg^ (1) 

where p is the proper energy density, p the pression and u M the unitary 4- velocity. The 
numerical current is 

jn = nuy. (2) 

where n is the proper density of particle. 

The motion of the fluid is governed by the equations of conservation 

V M T^ = (3) 

V,f = (4) 

The quantities - and £ are respectively the volume and the proper energy per particle. 
The enthalpy per particle is 

h= P -±V. (5) 

n 

The entropy s and the temperature T are defined in the thermodynamic relation 

dh = Tds H — dp (6) 
n 

so that h, T and - can be considered as functions of s and p. In particular h = h(s,p), n = 
n(s,p). It will be usefull to inverse the first relation and have 

p = p(s, h) , n = n(s, h) (7) 

Contracting the equation of motion (|3]) with u v gives 

V^nW*) -^8^ = 0. (8) 

Using (||), the equation (^) can be written 

nhu^V^ + « + u^u u )d^p = (9) 

Hence with (||) this equation is finally written 

8 h TB s 

u^,u u + « + u»u v ){^ £-) = (10) 
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and the equation (|8|) with (^j) gives 



ufdps = 0. (11) 



The entropy is conserved along a streamline. 

The following treatment is given by Lichnerowicz [0]. We suppose that the motion of 
the fluid is isentropic that is to say 

ds = (12) 

The equation @ is reduced to 

dh = —dp. (13) 
n 

In the relations ([?]) we can forget the dependence in s and have p = p(h), n = n(h). 
The equation flIU|) becomes 

u»V,u v + («# + it'X) (^) = (14) 
Let us introduce the quadrivector C 

C M = hii/j, (15) 
and its exterior differential Q = dC which is named vorticity. Its components are 

CV = V M a - V,C M (16) 

The equation (|I4]) can be written 



C% = (17) 

or zcfi = 0. 

The Lie derivative of Q along C 

C C VL = (di c + i c d)VL = (18) 

from the preceeding equation and dfl = 0. It is the Helmholtz theorem for an isentropic 
fluid: The vorticity is conserved in the motion. 

The motion is called irrotational if its vorticity is null Q = 0, that is dC = 0. If Q = 
on spacelike hypersurface, Q is null everywhere by the Helmholtz equation and the motion 
is irrotationnal. In this case the equations are 

V„C„ - V,C M = (19) 

V^f = (20) 



with 



J, = ^C, , h* = -C»C v . (21) 
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Let us note that the pression is p = j n{h)dh and that the velocity of sound v 
is given by 

1 n'(h)h 



[22] 



v 2 n(h) 

The invariance of the flow C M in the transformation generated by a Killing vector £ 

h 



is given by C^C = .This condition ensures also the invariance of since ^j^- is only 



dependent of h = ^J—C^C^ . If the fluid is irrotational we have 

C^C = (c% + i(.d)C = d(i^C) = (23) 

hence 

= constant (24) 
In a Minkowski space equiped with cylindrical coordinates 

ds 2 = -dt 2 + dz 2 + dr 2 + r 2 d(f) 2 , (25) 

let us consider a stationnary cylindrical solution of the equations ([19]) and (^). We call 
this motion a potential cylindrical vortex. The symmmetries are generated by three Killing 
vectors £ = J^, Jj, ^, so that 

C„ = (-£,L,0,M) (26) 

where E, L, M are the three contants of motion associated with the three vectors of Killing. 
In ( |2"E| ) we have also supposed that there is no source or sink in the core of the vortex. In 
consequence the radial component is null which ensures the equation (p0|) . There is yet 
an indetermination on n(h) which can be resolved by an equation of state. By a Lorentz 
transformation along z we can choose L = and have a circular flow. The scalar quantity 
h is calculated: 

M 2 

h 2 = -C^C^ = E 2 -L 2 --^. (27) 

We see from ( |2"T| ) that becomes spacelike in the neighbourhood of r = which is 
the central line of the vortex. So a central region, named the core of the vortex, must 
be empty which can be ensured by the equation of state of the fluid. We shall consider 
two different situations: The first is given by the familiar polytropic equation of state 0] 
p = kn 1 which can be rewritten using (|5]) and (|T3"| ) n = -^{h — m)~ . The quantity m is 

5 
3 



a constant of integration interpreted as the masse of a particle of the fluid. (7 is equal to 
or I following that the fluid is non relativistic or ultrarelativistic) . The second is the less 
usual stiff material where the the velocity of sound is equal to the velocity of light {v = 1). 
Integrating fl2"2"| ) gives n = kh. The radius r of the core is defined when the density n of 
the particules becomes null starting from the exterior. For the polytropic equation it is 
given by h(ro) = m and for the stiff material by h(ro) = 0. 
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3 Perturbation of the cylindrical vortex 



We consider now a little deformation of the relativistic cylindrical vortex described above. 
Its central line has an arbitrary shape but so that its curvature is weak and its deformation 
is slow. The section of its core stays circular with a little radius r . 

One introduces the local coordinates system (r A ,p a ) attached to the central line of the 
vortex: 

a? = X»{t a ) + p a N^(r A ) (28) 

The Minkowskian coordinates x M are expressed in function of the two coordinates r A = 
(t°, t 3 ) of the world sheet sweeped by the central line, and of two coordinates p a = (p 1 , p 2 ) 
pointing in a direction orthogonal to the world sheet along the two orthonormal vectors 
N£. We introduce also polar coordinates: 

p 1 = r cos <p , p 2 = r sin cf). (29) 

In this problem we need to define a reference length, choosed unity by convenience, in 
order to characterize the spacetime neighbourhood of a point of the central line in which 
we shall study the vortex. We shall suppose that the characteristic length I on which we 
have a change in the tangent plane of the world sheet is large and that the radius of the 
core r is small 

I > 1 , r < 1 

Let us note that we have two typical lenght I and r . 

The vetors N£ perpendicular to the world sheet and the tangent vectors |=pr have little 
change on a length unity. We can express this fact by 

r A dX^ r A 

K = K( T ) , -^a= x a( t ) (30) 

hence 

9aN^^ a , d B X^ = ±X%. (31) 
The Minkowkian metric g a p in the system of coordinates (r A ,p a ) can be expressed as 

9ab = Iab + 2K aAB p a + (k&Kobd + Wb) P a p b (32) 

9Ab = p a £ab^A (33) 
9ab = 5 ab (34) 
wherein we recognise the induced metric of the world sheet 

7AB = ^X>^ = 0(i), (35) 
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the extrinsic curvature 

KaAB = -^KaKb = 0{\) (36) 

and the twist defined by 

e ab u A = j V , u N£ A Nt; = 0(j). (37) 
Some other quantities will be useful: 



the mean curvature 



where 



and 



7 = det7AB, (38) 

K a = K aABl AB , (39) 
g = detg a p = ~/D 2 (40) 

D = 1 + K a p a + \{K a K h - K A B K* A )p a p\ (41) 

g AB = 1 AB_ 2K AB pa + 0( rJ ) (42) 

In the following we suppose that there is no twist in the deformation so that the cross 
terms of the metric cancel. 

It is natural to adopt the assumption that in the neighbourhood of the core whose the 
linear dimentions are some unities of r , the components C^, expressed in a orthonormal 
tetrad coincides with the solution of the stationary cylindrical vortex described in the 
preceding section, that is with (|26|) which in Minkowskian coordinates, is rewritten 

C U =(-E,L,-^M^ 2 M), (43) 

We can choose the coordinates r A on the world sheet so that jab takes the conformally 
flat form: 

Iab = F{—)r] AB (44) 

Let us note that OaF = O(j). The holonomic base 8a, d a is orthogonal but not orthonormal 
so that the identification of the solution with the stationary cylindrical solution must take 
into account the norm yF of 8a, hence 

C a = (-VFE,VFL,-^M,^M), (45) 

is the solution expressed in the holonomic base (8a, d a ). When expressed in the orthonor- 
mal tetrad (-j=,d a ) we can see that this solution coincides with the solution ( f43l) of the 
cylindrical vortex. 

We shall write the equations of motion fli~9|) and (|20D on the frontier of the core r = roand 
expand in power of the small parameter I -1 . The zero order which is the solution of 
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the straight vortex vanish identically. In this expansion appears also the second small 
parameter r . 

For this purpose it will be useful to rewrite the expression fl4"5|) in a more condensate 
form: 

C a = (C A VF,C a ), (46) 

where 

C A =(-E,L) , C a = {- P - 2 M, P - 2 M) (47) 

C a is of order rg 1 . 

We must calculate some quantities. Using fl4*2|), ([44]), the contravariant components C a 
are given by 

C A = g AB C B = -^=V AB C B - 2K AB p a VFC B + ... (48) 
V 

C a = 5 ab C b . (49) 



In (|48|), the second term is of order ^f-. The ellipse designates smaller terms; we shall follow 



this convention below. We can express 



M 2 

h 2 = -g aP C a C p = -g AB C A C B - — (50) 
in function of the corresponding quantity of the straight vortex which will be renamed 

hi = -V AB C A C B - (51) 

We obtain 

h 2 = h 2 + 2K AB p a C A C B F + ... (52) 

and 

W(h 2 ) = W(h 2 ) + 2K AB p a C A C B F^(hl) + .... (53) 

where we put W(h 2 ) = 

We can now expand the equations of motion on the frontier r = r$ of the core of the 
vortex. The equation (p0|) in coordinates (r A ,p a ) 

' d a (FDW(h 2 )C a ) = (54) 



FD 

gives 

± ^(%)ri™dWF + W(hZ)FK a C b fi* + ^(h 2 Q )F 2 K a C b 5 ab] j + . . . = (55) 

with 

K a = 2K c a D C c C D (56) 
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In this derivation we have used that for any function G(r): 



d a [G{r)5 ab C b ) = 

In the equation ([55]) the first term is of order j, the second and the third are in So 
we shall retain only the second and third terms in ^ which are the leading terms of the 
expantion (|55|): 

(w(hl)FK a + ^{hl)F 2 K^ C a = (57) 
It is easy to see that the leading term of the equation ([Tj]) is of order \. So the equation 



(|T9|) does not contribute at the same order as the equation (|20|) and can be forgotten. 

The equation fl5"7|) is written on the frontier r = r$ of the core; C a depends of the polar 
angle (p. Since the equation ( j57|) must be verified for all the polar angles when we turn 
around the core, we obtain finally the equation: 



W(h 2 )FK a + ^(h 2 )F 2 K a = 



(55 



Since the radius r of the core is supposed small, this equation gives a good approximation 
of the equation of motion of the central line. We shall discuss more extensively this equation 
in the next section. It is convenient now to examine the particular case of the stiff material 
where W{h 2 ) — f = k. We have |^ = and the equation flolf ) is reduced to 

K a = (59) 

which express the Nambu Goto dynamics of the central line in the intrinsic coordinates 

Let us project the extrinsic curvature K g ab in the Minkowski space: 

Kl B = -5 ab K aAB N£ = 8 A d B X° (<5£ - v P aV CD e p c e%) (60) 

whith 

4 - §. (61) 



From (|59|) and (|60|) the calculus of K^ B ^ AB gives the usual string equations of motion in 
Minkowski coordinates: 

' --dAiV^ll^X^) = (62) 



4 Discussion 

We cast the equations of motion (|58| ) in the interesting form: 



K aAB S AB = (63) 
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where 

S AB = 2^(hl)rj AG C cV BD C D + W(hl) v AB (64) 

For a weak deformation of the cylindrical symmetry (/ 3> 1), a rough approximation of 
(HI) S ives: 

e x c „ e v c ~ , e* ~ e z 3 ~ 1 (65) 
So, from fl5Q) the leading terms are tranverse (z = x, y): 

K\ B « d^X 4 , (66) 

and the equations of motion of the core are approximated by 

S AB d A d B X i = 0. (67) 



The equation fl22|) can be written 



" (68) 



dh 2 \v 2 J 2h 2 



If we assume that the fluid has a polytropic equation of state, then we have seen at the 
end of the second section that when r tends towards the radius of the core r , then h, resp. 

dW 
dh 2 



n(h), tends towards m, resp. zero. As a result, the quantities W and 4p- tends towards 



zero but the quotient of ^ by W tends towards infinity. Therefore in the vicinity of r = tq 
we can neglect the second term in ( |B"4"D and after simplification the equation Q8TD is written 

E 2 8lX l + 2ELd d 3 X i + L 2 djX l = (69) 

This equation is parabolic and cannot represent propagating waves along the core of the 
vortex. To be more specific, we have seen in the second section that it is possible by a 
Lorentz transformation along z, to choose L = 0. In this case the equation ( |69D becomes 
particularly simple. Its solution 

X = a(r 3 )r° + 6(r 3 ) 

moves off the equilibrium linearly with the time. It can represent the motion of the vortex 
only at the beginning during a short time, since the method of resolution supposes a small 
deviation from the initial equilibrium, which is not the case at later time. 

The equation ( |62D can be expressed in the preceeding approximation to compare with 

^ d 2 X l - d$X { = 0. (70) 

A small perturbation can propagate along z. This result singles out the dynamics of Nambu 
Goto obtained in stiff material. 
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5 conclusion 



If we refer to the usual experience in hydrodynamics, it can appear somewhat strange 
that, in a polytropic fluid, a small deformation of the cylindridal core of the vortex cannot 
propagate as along an elastic string. May be the model of isentropic irrotational fluid is too 
restictive. We must also remember that in non relativistic fluid, if we bend a cylindrical 
line vortex, we induce infinite quantities which are dynamically embarassing [ |T0|j . 

On the other hand in stiff material the core of the vortex has the dynamics of the 
Nambu Goto string and consequently small deformations can propagate. 

To conclude, at our knowledge, it is the first time that the Nambu Goto dynamics is 
directly connected to the classical relativistic hydrodynamics. 
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